Interfaces of polydisperse fluids : surface tension and adsorption properties 
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We consider a system of spherical colloidal particles with a size polydispersity and use a simple van 
der Waals description in order to study the combined effect of both the polydispersity and the spatial 
non-uniformity induced by a planar interface between a low-density fluid phase (enriched in small 
particles) and a high-density fluid phase (enriched in large particles). We find a strong adsorption 
of small particles at the interface, the latter being broadened with respect to the monodisperse 
case. We also find that the surface tension of the polydisperse system results from a competition 
between the tendancy of the polydispersity to lower the surface tension and its tendancy to raise the 
critical-point temperature (i.e. its tendancy to favor phase separation) with the former tendancy 
winning at low temperatures and the latter at the higher temperatures. 



I. INTRODUCTION 



Many of the complex fluids used in the industry or 
in the soft condensed matter physics laboratory are col- 
lections of nearly identical particles which exhibit one 
or several polydispersities, i.e. properties such as the 
size or the shape of these particles which are distributed 
whithin some interval in an almost continuous manner 
[1]. These fluids are hence continuous mixtures of sim- 
ilar particles and it is of some practical importance to 
know how their composition or polydispersity influences 
their physical properties, e.g. their phase behavior and 
rheological properties. Here we will be concerned only 
with the equilibrium phase behavior of such polydisperse 
fluids. The generalization of the well-established meth- 
ods for the study of phase transitions in discrete mixtures 
to continuous mixtures is a technically very demanding 
task which has recently become an active area of research 
[2]. Most of this research has been limited to spatially 
uniform (or bulk) phases while it is our purpose here to 
extend it further to non-uniform situations involving the 
interface between two coexisting bulk phases. Such an 
intrinsic interface is different from, e.g. the interface be- 
tween a bulk phase and a substrate (see [3] for some pre- 
liminary work in this direction). Indeed, in the latter case 
the interface can be characterized thermodynamically by 
a surface excess free- energy defined relative to the sub- 
strate whereas in the former case it is characterized by a 



surface tension (which is both a surface free-energy and a 
thermodynamic force) defined with respect to an intrin- 
sic "surface of tension" [4] . We will hence be particularly 
interested in the influence of the polydispersity on this 
interfacial surface tension. 

In order to keep the problem manageable we will re- 
strict ourselves to the fluid phases of a system of spher- 
ical particles with a size-distribution, a situation typical 
for many colloidal dispersions [5]. The initial or parent- 
phase size-distribution will be assumed to be fixed once 
and for all by the production process of the colloidal 
particles and taken to be of the monomodal type, i.e. 
peaked around a single reference species, such as is ap- 
propriate for the polydisperse generalization of a one- 
component system. When this (initial) parent-phase is 
put into appropriate thermodynamic conditions it will 
phase separate or " fractionate" into two (or more [6]) 
daughter-phases. Since these daughter-phases differ in 
density and composition an interface will build up be- 
tween them across which the properties of one bulk phase 
transform continuously into those of the coexisting bulk 
phase. The properties of such a spatially non-uniform 
two-phase system are most conveniently studied in two 
steps. First, one determines the two spatially uniform 
bulk phases which are able to coexist in equilibrium. 
Next, one determines the profiles across the interface of 
those properties which are spatially varying in the two- 
phase system. For the first step we will use the results of 
our earlier study [7] based on the van der Waals (vdW) 
free-energy of a polydisperse system of spherical colloidal 
particles interacting via excluded volume repulsions and 
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vdW-like attractions. Such a description is of course not 
exact but simplifies considerably the technical problems 
raised by the study of phase equilibria in polydisperse 
systems. For the second step we will use an earlier ex- 
tension of the vdW free-energy to spatially non-uniform 
systems [8] and generalize it here to polydisperse systems. 
In order to extract the surface tension we will, finally, 
adapt to the present case a general procedure advocated 
elsewhere [9]. 

In Sec. H we introduce the vdW density functional of a 
spatially non-uniform polydisperse system. The density 
profiles across a planar interface are computed in Sec. Ill 
while the corresponding adsorption properties will be dis- 
cussed in Sec. [V. The pressure profiles across the planar 
interface are determined in Sec. ^ while the resulting 
surface tension is presented in Sec. Vl. Our conclusions 
follow in the final Sec. VII. 



II. THE SPATIALLY NON-UNIFORM 
POLYDISPERSE SYSTEM 

The equilibrium properties of spatially non-uniform 
systems (e.g. interfaces) are most easily studied within 
density functional theory (DFT) [10]. The starting point 
of DFT is the variational free-energy, A{T, [p], [<&]) : 



FiT, [p]) 
J dr J da p{r, a){^{r, a) 



K^)} (1) 



where T is the equilibrium temperature, F{T, [p]) the 
intrinsic Helmholtz free-energy viewed as a functional 
(indicated as [p]) of the average local number density, 
p(r, a) , r being the position variable (assuming spher- 
ical particles) and a the (dimensionless) polydispersity 
variable, fJ.{a) is the chemical potential of species a and 
$(r, cr) the one-body external field responsible for the 
spatial non-uniformity of the system (the functional de- 
pendence of A on $(r, cr) being indicated as [<&]). The 
equilibrium density, p(r,a), corresponding to a given 
$(r, cr), can then be obtained by solving the Eulcr- 
Lagrange equation : 



SA{T, [p], [<&]) , 
6p{r,a) 



= 



corresponding to (fil), viz. 



Ai(cr) = $(r,cr) 



6p{r,a) ' 



(2) 



(3) 



Eq. (|^) expresses the fact that in equilibrium the chem- 
ical potential fJ.{a) of each species a has to remain con- 
stant in space. In the present study, F{T, [p]), will be ap- 
proximated by the following vdW-type expression [7,8,11] 



^^(r, [p])^kBT J dr J dapir,a) 
A^{a)p{r,a) 

+ ^ J dr J da J dr' J da' p{r,a) 
VA{\r~r'\;a,a')p{r',a') 



(4) 



where ks is Boltzmann's constant, A{a) the thermal de 
Broglie wavelength of species a, VA{r; a, a') the poten- 
tial of attraction between two particles of species a and 
a' a distance r = \r\ apart, while E{r, [p]) represents the 
excluded volume correction resulting from the repulsions 



^(r, [p]) = 1- J dav{a)p{r,a) 



(5) 



v{a) = ^R^{a) being the volume of a spherical particle 
of radius R{a). It should be noted that, as is usual in this 
context [11], a is used here both as a species label and as 
the (dimensionless) polydispersity variable, R{a)/R{l), 
R{1) being the radius of the reference species a = 1. 
From (^) we obtain for (^ : 

M.) = ^(r,a) + fc.rin( ^yf;;) ) 



v{a) 



Eir, [p]) 

i.(r,[p])/^^'^(^''^') 
jdr' J da'VAi\r~r'\;a,a')pir',a'). 



(6) 



The above represents a straightforward extension of DFT 
to continuous mixtures while (@-||) reduces to the vdW 
free-energy used in [7] for the uniform polydisperse sys- 
tem as well as to the vdW free-energy used in [8] for the 
non-uniform monodisperse system. Both approximations 
reducing to the usual vdW free-energy for the uniform 
monodisperse system. The basic physics of the vdW ap- 
proximation being, as usual, the correction of the ideal 
gaz behavior for the finite size of the particles via the 
excluded volume correction (E) and the inclusion of the 
cohesion between the particles via the interparticle at- 
tractions (Va), as described here respectively by the first 
and second term of (^. Of course, more involved expres- 
sions of F{T, [p]) are available but these can only add 
further complications to the already fairly complex cal- 
culations required by the present combination of the non- 
uniformity with the polydispersity of the system. Past 
experience has shown however that the present vdW ap- 
proximation is able to capture the essence of the underly- 
ing phase behavior in a qualitatively correct manner [12]. 
In [7] we have studied several model-polydispersities dif- 
fering in the cr-dependence of v{a) and VA(r;a,a'). It 
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was found there that the model based on the simple ap- 
proximation : 

v{a)^v{l), VA{r;a,a')= a a' VA{r- 1,1) (7) 

has a phase behavior which is similar to that of the more 
involved expressions but is simpler to study. Hence- 
forth we will use thus (|-|) together with (@). The 
physical contents of reflects the fact (cf. [7]) that 
the amplitude-polydispersity of VA{r;a,a') dominates 
the volume-polydispersity of v{cr). The inclusion of the 
volume-polydispersity of v{cr) will therefore not alter 
qualitatively our conclusions. 



III. THE PLANAR INTERFACE 

We will consider a planar interface perpendicular to the 
z-axis. Translational invariance in the (x,y)-directions 
implies then, p{r,a) p{z,a) and $(r, cr) — > $(z,cr), so 
that the Euler-Lagrange eq. (||) can be rewritten, after 
separating the local (in z) and non-local contributions, 
as : 

^((t) = cr) + po{z, cr; T, [77]) 
+a dz'Vi{\z-z'\){m{z')-miz)} (8) 



with po{z, a; T, [rj]) a shorthand notation for : 

lio{z,cr;T, [rj]) = fcg T In ^ + fcg Tin 

v{l) l-?7o(z) 

+fcsr-^^^+aT/om(z) (9) 
1 - Vo[z) 

where 

r]{z,a) ^ v{l) p{z,a), r]o{z) ^ J daij{z,a), 

771(2;) = / d<7arj{z,a), (10) 



are the dimensionless density and polydispersity mo- 
ments, whereas : 



00 />oo 



v{l)Vi{\z\)^ / dx / dyVA{r;l,l), 

J —oc J —00 

vil)Vo= I dry^(r;l,l), (11) 



with (7=1 denoting the reference particle of volume 
v{\). The external (symmetry breaking) fleld, $(z,cr), 
will as usual be replaced by boundary conditions. We 
thus consider eq. (H) without external field, viz. : 



p[a) ^ PQ{z,(j;T,[ri]) + a / dz' Vi{\z ~ z' \) 

JO 

{^,{z')-r„{z)} (12) 



and require that for, z ±00, the solution r}{z,a) of 
( |T^ ) matches the bulk- phase densities, say 7y±(cr), or 
'q{z = ±00, cr) = rj±{a). These bulk-phase densities must 
hence satisfy eq. ( p^ for z — ±00. Taking the limit of 
( p^ ) for z — > ±cxD, the second term in its r.h.s. will vanish 
and we obtain : 



p{a) = po{z = ±co, cr; T, [r/±]) 



(13) 



i.e., the chemical potential of species cr in the non- uniform 
system must be constant and equal to the chemical po- 
tential of species a in the two bulk-phases. Indeed, eval- 
uating the r.h.s. of ( p^ ) from (|^) for 77(2 = ±00, cr) = 
77± (cr) we obtain : 

MO (±00, (i; T, [,y±]) = fcs T In + fcs T In ^'^^ 



v{l) 



4 



+kBT-^ + c7Vovf (14) 
1 



where : 



(15) 



while the r.h.s. of ( p^ ) represents (cf. [7]) the chemical 
potential of a uniform phase of density 77^ (cr) . When the 
two bulk phases are in equilibrium, the chemical poten- 
tial of the 77+ (cr) phase must be equal to that of the 77- (cr) 
phase, hence ( |l3|) will be satisfied. Eq. ( |l3|) allows us to 
elimite /i(cr) from (|l^) and rewrite it as : 



Mo(±oo, cr; T, [?7±]) - fiQ{z, a; T, [rj]) 
dz'V^{\z-z'\){m{z')-miz)} 



(16) 



an integral equation for 77(2, cr) incorporating the bound- 
ary conditions. On using (^ and ( |l4|) we can rewrite ([l^ ) 
as : 



77(z,a)=A±(z)Af±(z,a) 
where A^(z) is a shorthand notation for : 

A-iz)^'-^eM ' 



1 



4 



1 



1 - Vo{z) 



(17) 



} (18) 



and M^{z,a) for : 
Af^(z, cr) — ri±{a) exp a 



dz'/3yi(|z-z'|){77±-77i(z')} 
(19) 



where /? = l/ksT. Taking now the first two cr-moments 
of (|r^ yields a system of two integral equations for r/o(z) 
and ?7i(z), viz. : 
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r^^{z) = A^{z)M^{z) 
77i(z) = A±(z)M±(z) 



(20) 



where 



M^{z)^ J dcrM^{z,a), M^{z) ^ J daaM^{z,a). 

(21) 

Solving (|o|) and substituting the result into ( p^[jT9| ) 
yields finally ri{z^a) via (p^. Note that, since the bulk- 
phase densities rj± (cr) must correspond to the same chem- 
ical potential (cf. (Esl)) : 



Mo(oo,cr;T, [r/+]) = ^o(-oo,(t;T, [ry_]), 



(22) 



the equations (p^pO|) bearing the sign are equivalent 
(but not identical) to those bearing the "-" sign. The 
above equations can therefore be used with cither sign, 
the results will be the same provided, of course, that the 
two bulk phases are coexisting equilibrium phases. 

To proceed we must specify VA(r; 1, 1). Since only fluid 
phases are involved the particular form of Va{t", 1, 1) is 
not very important and in view of (|ll| ) we will, for sim- 
plicity, take it to be gaussian or in the notation of [7] 



VA{T;l,l) = -e{l,l)%v{\) 



exp(- 



^3(1) (1)3/2 



R{1) 
(23) 



where e(l,l) is a reference amplitude. If, t — 
ksT / t{l,l), denotes the dimensionless temperature eq. 
(O) yields on using (P3|) : 



(iVi{\z\) 



8 exp(— 6z^) 



R{1) 



t 

(24) 



Below we have used (for convergence reasons) h = A. 

Finally, eq. ( ^7|) also requires explicit data for ri±{a). 
For the latter we will take the two-phase coexistence den- 
sities obtained in [7] for the same temperature {t) and for 
an initial parent-phase density, pQ{cr) = pohQ{a), of av- 
erage density po (or, in dimensionless form, rjo ~ pQ v{l)) 
and a Schulz-Zimm size-distribution ho{a) : 



'^exp[—aa) 



(25) 



1 



where T{a) is the Euler gamma function and / 
(1/a) the polydispersity index. Note that 1/a = / — 1 
is the variance of ho{a) so that 1 = 1 (or a = oo) corre- 
sponds to the monodisperse limit whereas the reference 
species a = 1 corresponds to the average value of a in 
the parent-phase. Of course, other size-distributions can 



be used but as shown in [6,7] the particular form of ho{a) 
has little influence as long as it remains monomodal. 

In the present work, we have studied two polydisper- 
sities, viz. a = 50 (/ = 1.02) and a = 15 (/ = 1.07), for 
several temperatures {t) and densities (770). In Fig. 1 we 
show three binodals of the bulk phase diagram for a = 50 
(see also [7]). In Fig. 2 we show the size-distributions 
of two bulk-phases coexisting for t = 1, rjo = 0.48 and 
a = 15. The inset of Fig. 2 shows the corresponding den- 
sity distributions. It is seen there that for a range of a- 
values (0.3 < a < 0.7) the densities of the "low"-density 
phase actually exceed the corresponding densities of the 
"high" -density phase. Finally, in Fig. 3 we show a variety 
of density profiles for the non-uniform two-phase system. 
It is seen that, compared to the monodisperse case, the 
polydispersity widens the interfacial region. The results 
of Fig. 3 have been obtained by solving eq. (|2^) it- 
eratively (e.g. by starting from a tanh-proflle) whereas 
the results shown in Figs. 1-2 have been obtained as 
explained in [7]. 



IV. ADSORPTION PROPERTIES 

In macroscopic thermodynamics it is customary to 
replace the continuous density proflles, ri{z,a) = 
w(1)p(z,(t), obtained in Sec. ^ by discontinuous (with 
respect to z) profiles, r]{z,a) = v{l)'p(z, a), of the form 
[4] : 

p{z, cr) = P+{cr) 9{z - zg{(j)) + P-{cr) 0{zg{(j) - z) 

+V{a)5{z-ZG{o)) (26) 

where ?7±(cr) = t;(l) p±{<t) are the density distributions of 
the two coexisting bulk phases and V{(j) is the adsorption 
of species a at the interface for which z = ZQ^a) is the 
Gibbs dividing surface of species cr. In (p^), 9{z) denotes 
the Heaviside step function and 6{z) the Dirac delta func- 
tion. The macroscopic (^(z, cr)) and microscopic (p(z, cr)) 
profiles can be adjusted by requiring them to satisfy : 



dz {p{z, cr) — p(z, cr)} = 



(27) 



which implies that r(cr) be defined as the surface excess 
density, viz. : 



Tia)= / dz{p{z,a)-p{z,a)} (28) 

J —oo 

where p{z, cr) is the following bulk-phase switch function 

p{z, cr) = p+{a) 9{z - zg{(j)) + p-{cr) 0{zg{o-) - z). 

(29) 

Since at z = ±00 both densities match, p{z = ±cx), cr) = 
p{z = ±cx),cr) = p±{<y), we can integrate (28) by parts 
and obtain : 
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r(a-) = / dz {zg{(t) ~ z) p'{z,a) 



(30) 



where p'(z,cr) = dp{z,a)/dz. As seen from ( |30| ) the 
value of r(cr) attributed to a given p{z, a) stih depends 
on the value of zg(o'). Since there is no absolute deter- 
mination possible for zcicr) we have to fix it arbitrarily, 
e.g. for the reference species a = 1. Taking henceforth 
zg((t) = zg{1) for all a we can fix 2g(1) by requiring 
that the corresponding adsorption, r(l), vanishes. Eq. 
(^) implies then : 



^g(1) = 



J^^dzzp'jz,!) 



whereas (pOh becomes : 



ri(cr) = / dz{zG{l) - z) p'{z,a) 



(31) 



(32) 



where the subscript 1 on T(a) indicates that the adsorp- 
tion of species a is referred to the zero-adsorption Gibbs 
dividing surface ( |3l| ) of the reference species cr = 1, hence 
ri(l) = 0. Since, moreover, the system is of infinite ex- 
tend in the z-direction we may choose zg(1) as the origin 
of our coordinate system, i.e. zg(1) = 0. Some examples 
of ri(CT) are given in Fig. 4. As seen from Fig. 4, at 
the interface there is both an excess of small particles 
(Fi (cr) > for cr < 1) and a depletion of large particles 
(Fi(cr) < for cr > 1) with an adsorption (Fi(cr)) which 
strongly depends on t and r/g. 



V. PRESSURE PROFILE ACROSS THE PLANAR 
INTERFACE 

Besides the density profile {p(z, cr))which gives rise to 
the adsorption properties described in the previous sec- 
tion, an interface also involves a pressure profile {p{z)) 
which in turn gives rise to the surface tension as will 
be shown [9] in the next section. In order to expose 
the pressure in the interior of the interface described by 
p{z,a), we first cut this interface with a plane perpen- 
dicular to the density profiles, say the x = plane, and 
remove the matter on the x < side of this plane while 
leaving the matter on the x > side intact. Such a semi- 
infinite system with a density, p{r, cr) ~ 0{x) p{z, a), can 
be realized within the DFT of Sec.|| by replacing the 
matter removed from the a: < half-space by a corre- 
sponding external field, say <i>(r, cr). The pressure act- 
ing normal to the a; = plane, i.e. acting in a di- 
rection which is tangential to the density profiles, can 
then be obtained by submitting the x = plane to 
an infinitesimal non-uniform normal deformation, viz. 
X ^ X + 5u{y, z), and computing the resulting thermo- 
dynamic work of deformation (cf. [9]). Since during this 
infinitesimal deformation, $(r,cr) ^{r,a) -|-(5$(r, cr) 



and p{Y,a) — > p{r,a) + Sp{r,a), the system has to re- 
main in equilibrium at the given T and p{(t), the relation 
between 5^{r, a) and Sp{r, a) can be obtained from the 
equilibrium condition (tBI) as : 



(5/x(cr) = = (5$(r, cr) + y dr' J da' 

— , , JT Op{r , cr ). 

6p{r,a) dp{r',a') 



(33) 



The resulting thermodynamic work of deformation, SA, 
can then be obtained from (|l|) : 

^ JdrJ dap{r,a)S^r,a) (34) 

or on using (|33|), from : 

^^|t,[^] = — J dr J da J dr' J da' p{r,a) 

— — — — -5p[r,a). (35) 

5p{r,a) Spir'^a') 

Since in the present geometry we have, 5p{r, a) — p(x + 
Su{y, z), y, z, a) - p{x, y, z, a) = 6u{y, z) .{dp{r, a)/dx) + 
0{6u^), eq. ( |3^ ) can be rewritten after dropping the 
0{6u^) term: 



dy' / dz' 5u{y' ,z')p{y' ,z') 



(36) 



which defines the pressure p(t/, z) acting at r = (0, y, z) 
in a direction normal to the x = Q plane. Indeed, since 
5u{y^ z) is arbitrary (p5mq) imply : 



(37) 



p{y,z) — J dx J da J dr' J da' p{r' ,a') 
6^FiT,[p]) dp{r,a) 



6p{r' ,a') Sp{r,a) dx 



where, for convenience, we have interchanged the role of 
the primed and unprimed variables. In the present vdW- 
approximation we obtain from ^ : 



S'FjT, [p]) 
Sp{r,a) 6p{r',a') 



6{r-r')6{a~a') 
P(r,cr) 

VA{\r~r'\;a, a') 
Sir - r') 



Eir, [p]) 

which on behalf of (|^), reduces here to 



S'FjT, [p]) 
Sp{r,a) Sp{r',a') 



knT 



dir^r')S{a-a') 
P(r,cr) 
aa'VA{\r~r'\-l,l) 

-f^BTSir-r'Ml)ii±^} 
{E{r, [p])y 



(39) 
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so that (pf\j can be rewritten : 



p{y,z)^ j dx J daiksT 
J dr' J da' p{r',a')VAi\r-r'\;a, a 



d p{r,a) 
dx 



dx 

+kBTv[l) p{r,a)—— — / da }. (40) 



dx 



Taking into account that here, p{r,a) = 9{x) p{z,a), we 
can rewrite ( ^0| ) as : 

J-oo clx l-rio{z)6[x) 
+ f dx f dr'VA{\r-r'\;l,l)S{x),^i{z)e{x')r]i{z') (41) 



where ?7o(z) and rji^z) have been defined in (10). Eq. 
( ^ ) can be rewritten as : 

v{l)p(z) = t^lM^ + I ^,(,) r dz'V,{\z - z'\) vAz') 

(42) 

where T^d^^l) was defined in ( pi] ) and we took into ac- 
count that p{y, z) is independant of y as expected from 
the translational invariance in the y-direction. We finaUy 
rewrite ( p2[ ) in a manner similar to (|8|) : 

v{l)p{z) = ■y(l)po(2;7', [?7]) 

1 /"^ 

+ -r;i(z) d^Vi(|z-z'|){r7i(z')-r7iW} (43) 
with, po{z;T, [rj]), a shorthand notation for : 



v{l)po{z;T, 



'-f^ + lvoimizyr (44) 

l-?/o(z) 2 



where Vq was defined in (|T^). It is seen that (^) repre- 
sents the usual vdW-pressure of a uniform (polydisperse) 
system evaluated for the local density ri{z,a), while the 
same is true of (||) for the chemical potential. From ( ^ ) 
it is seen that the local pressure p{z) is completely de- 
termined by the local density ry(z, a) . From ( ^3| ) and 
ri{zLoo,a) = ri±{a) we obtain p{±oo) — pq{±oo;T, [r]±]), 
but since rj±{a) must satisfy [7] : 



kBT% . 1 



ksTriQ 1 



+ 77 K, [vtY = -i^^ + -,V, iv^y (45) 



1 - Vo 2 



together with (p^), we have p(oo) = p(— oo), or 
/ dzp'{z) — 0, which expresses the stability of the planar 
interface (cf. [9] for details). 

Some of the pressure profiles, p(z), obtained from (^) 
using the density profiles of Sec. Ill are shown in Fig. 
5. While p{z) remains constant in the bulk phases it ex- 
hibits a structure in the interfacial region which is more 



pronounced for the lower temperatures and disappears 
gradually when t = tc '^s approached. This structure con- 
sists of a pressure depletion on the high-density side of 
the interface and a pressure excess on the low-density 
side (a similar structure was found in the monodisperse 
case [8] for a different potential Va(?'; 1, 1) and was seen 
also in the simulation results of [13]). This local structure 
of p{z) reflects a competition between the characteristic 
length-scales of p{z, a) and of VA(r; 1, 1). It is also seen 
(compare Figs. 5(a) and (b)) that increasing the poly- 
dispersity (i.e. lowering a) widens the interfacial region 
(i.e. the region where p'{z) ^ 0). 



VI. SURFACE TENSION AND SURFACE OF 
TENSION 



In a way analoguous to (p6[), the microscopic pressure 
profile p{z) of Sec. ^ is replaced in macroscopic thermo- 
dynamics by a discontinuous pressure profile, p{z) : 

p{z) = p+ 9{z - zo) + p_ 0{zq - z) - -1 5{z - zo) (46) 

where p± — p(±oo) denote the bulk-phase pressures and 
7 is the surface tension acting on a surface of tension 
located aX z — zq. As in (|2^) the two profiles can be 
adjusted by requiring that: 



dz {p{z) — p{z)} = 
which on behalf of (|4^) yields : 

/oo 
dz {p{z) -p{z)} 
-oo 

where p{z) : 



Zo) 



P- f^[ZQ 



(47) 



(48) 



(49) 



is the switch function for the bulk-phase pressure (cf. 
(p9|)). Integrating ( ^ by parts and taking into ac- 
count that for a planar interface we must have, p{ioo) — 
p{±oo) = p± together with — p_, yields: 



7 = 



dz zp'{z) 



(50) 



which shows that 7 can be determined from the knowl- 
edge of p' {z) — dp{z)/dz alone, i.e. without knowing zq 
(cf. the difference with (|30|)), a feature specific to the 
planar interface (cf. [9]). To determine zq one can nev- 
ertheless also impose (cf. [9]) that : 



(51) 



dz z {p{z) — p{z)} = 
which on using (H^) can be rewritten as : 
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= 



dz [z - zq) {p{z) -p{z)} 



dz [z - zo) {p{z) ~p{z)} 



(52) 



or integrating (52) by parts and using — p_ one ob- 
tains : 



dz {z — zo)'^ p'{z) = 0. 



Using, / dzp'{z) = 0, eq. (p3[) yields finally 



^ 1 I^oodzz^p'jz) 
2 jT^dzzp'iz) 



(53) 



(54) 



so that both 7 and zq can be obtained from p'iz) (of. 
(1^) and (H)). 

Fig. 6 shows some of the results obtained for 7 us- 
ing the pressure profiles of Sec. |^. From Fig. 6a it is 
seen that 7 decreases when increasing 770 at a fixed poly- 
dispersity. When, instead, the polydispersity is changed 
starting from the monodisperse case {a — 00) there are 
two competing effects (cf. Fig. 6b) : the polydispersity 
tends to lower the surface tension but at the same time it 
raises the critical-point temperature since the polydisper- 
sity favors the phase separation (cf. [7]). As a net result 
of this competition the surface tension of the polydis- 
perse system is lower than that of the monodisperse one 
for t < tQ but exceeds it for t > Iq with a crossover tem- 
perature to which increases with the polydispersity. To 
show this more clearly Fig. 6c displays 7 versus t/tc{a), 
where tc{a) is the critical point temperature of the sys- 
tem with polydispersity a. 

The above surface tension (7) acts on the surface of 
tension located ai z = zq (cf. eq. (^6|)). From Fig. 7 
it is seen that zo < 0, i.e. the surface of tension is dif- 
ferent from the zero-adsorption Gibbs dividing surface 
{z — zg(1) = 0) and located on the high-density side of 
this interface (this can also be seen from Fig. 5). The 
fact that Zq 7^ zg{1) points to a fundamental inadequacy 
of the macroscopic description of interfaces (cf. [9]) since 
two different quantities are used to locate the interface 
on a macroscopic level. The quantity, It = zci^) — zq, is 
usually called Tolman's length [4]. Finally, from Fig. 8 it 
is seen that on approaching the critical point 7 vanishes 
with a classical critical exponent (= 3/2) as expected 
from the mean-field vdW-theory [4] . 



VII. CONCLUSIONS 

We have studied the planar interface resulting from 
the phase separation or fractionation of a parent phase 
of a polydisperse colloidal system into a low-density fiuid 
phase enriched in small particles and a high-density fluid 



phase enriched in large particles. In order to tackle the 
combined effect of the system's polydispersity and spatial 
non-uniformity we have kept the theorical description as 
simple as possible but are confident that similar results 
can be found from more involved descriptions. Using a 
simple van der Waals description [7-8] to model the poly- 
disperse non-uniform system of spherical colloidal parti- 
cles with excluded volume repulsions and gaussian at- 
tractions it was found that the small particles accumu- 
late at the interface, the latter being moreover depleted 
with larger particles and broadened with respect to the 
monodisperse case. We also found that for a given tem- 
perature the surface tension is the result of a competi- 
tion between two polydispersity-induced effects, namely 
its tendancy to lower the surface tension and at the same 
time to raise the critical point temperature, with the for- 
mer effect winning at low-temperatures and the latter 
at higher temperatures. Finally, the surface of tension 
was found to be located on the high-density side of the 
interface pointing to a positive Tolman's length. 
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Figure Captions 



FIG. 1. The temperature (t)-density (77) bulk-phase di- 
agram for a = 50 (cf. [7]). Three binodals are shown. 
They correspond to the parent-phase densities i]q ~ 0.3 
(dashed hne), ?7o = '7c = 0.3659 (full line) and ?7o = 0.45 
(dot-dashed line). Two binodals are truncated upwards 
at resp. the supra-critical temperature t ~ 1.246 (for 
r?o = 0.3) and the infra-critical temperature t ~ 1.1905 
(for 770 = 0.45) while the untnmcated binodal passes 
through the critical point tc — 1.2355, r^c = 0.3659. Sim- 
ilar results are obtained (cf. [7]) for a = 15 in which case 
the critical point corresponds to = 1.2889, 77c = 0.4842. 

FIG. 2. The density- (r;± (a) = r]±h±_{o)\ cf. inset) 
and size-distributions (/i±(cr)) of the low-density (77+; full 
line) and the high-density (77-; dashed line) bulk phases 
(77+ < 77_) which coexist for i = 1, 770 = 77c = 0.4842 and 
a = 15. Note from the inset that for, 0.3 cr < 0.7, the 
behaviors of rj±{a) and r]± = J dar]±{a) are reversed, i.e. 
although 77+ < rj- we have 77+ (a) > ri-{a) for these a- 
values. The low (high)-density phase is enriched in small 
(large) particles. 

FIG. 3. Density profiles for a = 15 across a pla- 
nar interface : (a)-(b) at constant temperature (t — 1) 
and (c) at constant density {rjo = ijc = 0.4842). In (a) 
we show 77o(^) (full hne) and 771(2) (dotted line). Also 
shown for comparison is the monodisperse case (dashed 
line) corresponding to a = 00 and 770(2;) = 771(2). It is 
seen that the polydispersity broadens the interface. In 
(b) we show rj{z,a) for a = 1.25 (dot-dashed), tr = 1 
(dashed), cr = 0.75 (dotted) and a = 0.65 (full line). It 
is seen that the small particles (cr < 1) accumulate in the 
interfacial region. Note also that 77+ (0.65) > 77_(0.65) 
whereas 77+ (0.75) < 7?_(0.75) in agreement with the re- 
versal seen in Fig. 2. In (c) we show ri{z, a = 0.65) for 
a = 15, 7^0 = 0.4842 and t = 0.85 (full line), 1 (dots), 
1.1 (short dashes), 1.2 (long dashes) and 1.28 (dot-dash). 



(Here z* = (z - zg(1))M1)). 

FIG. 4. The adsorption ri(cr) of cr-particles relative 
to the zero- adsorption Gibbs dividing surface (zai'i-)) of 
the reference particle (a = 1) at (a) fixed density and 
(b) fixed temperature. Panel (a) corresponds to a = 15, 
,70 = 0.4842 and t = 0.90 (full line) 1.00 (dots), 1.10 
(short dashes), 1.20 (long dashes), 1.28 (dot-dash). Panel 
(b) corresponds to a = 50, f = 1.15 and 770 = 0.45 (hne), 
0.3659 (dots), 0.3 (dashes). Note the rapid variations 
with t and rjQ of the interfacial excess of the small parti- 
cles {a < 1) and depletion of the large particles {a > 1). 
(Here rt(a) = ri(a). v(l)/ii(l).) 

FIG. 5. The pressure profile (p(z)) versus the distance 
{z) from a planar interface perpendicular to z-axis for (a) 
a = 50 and 770 = 0.3659 and (b) a = 15 and 770 = 0.4842 
and three temperatures t = 0.9 (full line), t = 1 (dots) 
and t = 1.1 (dot-dashes). The interfacial region is seen to 
be broadened by the polydispersity. All profiles exhibit 
a pressure depletion (excess) on the high (low) density 
side of the interface. (Here p*{z) = p{z) i'(l)/e(l, 1) and 
z* = z/R{l).) 

FIG. 6. The surface tension (7) versus the tem- 
perature {t) for : (a) a = 50 and 770 = 0.3 (dots), 
110 = 110 = 0.3659 (full line) and r]o = 0.45 (dot-dash); 
(b) a = 15 (dot-dash), a = 50 (dots) and a = 00 (full 
line) for 770 = ric{a); (c) the same as (b) but plotted now 
versus t/tc{ct) where tc{cx) and r]c{a) are, respectively, 
the (reduced) critical-point temperature and density of a 
system with polydispersity index 1=1 + (1/a). (Here 

7*=7Kl)Ml,l))i?(l)-) 

FIG. 7. Tolman's length {It = zg(1) — zq) versus 
the temperature (t) for a = 50 and tjo = 0.3 (dots), 
110 = 110 = 0.3659 (dash) and 770 = 0.45 (dot-dash). (Here 

= It/R{1)-) 

FIG. 8. A to 7* versus Int* plot for a = 00 (fuU 
line) and a = 15, rjo = Vo = 0.4842 (dot-dash). In both 
cases the critical exponent (=3/2) of 7 is classical. (Here 



9 



Fig. 1 , Bellier-Castella et al, PRE 
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